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ABSTRACT 

The Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory (SDO) has observed 
two sun-grazing comets as they passed through the solar atmosphere. Both passages resulted in a 
measurable enhancement of Extreme Ultraviolet (EUV) radiance in several of the AIA bandpasses. 
We explain this EUV emission by considering the evolution of the cometary atmosphere as it interacts 
with the ambient solar atmosphere. Molecules in the comet rapidly sublimate as it approaches the 
Sun. They are then photodissociated by the solar radiation field to create atomic species. Subsequent 
ionization of these atoms produces a higher abundance of ions than normally present in the corona 
and results in EUV emission in the wavelength ranges of the AIA telescope passbands. 
Subject headings: Comets: general — Comets: individual (C/2011 N3) — Comets: individual (C/2011 
W3) — Sun: corona — Sun: general 



1. INTRODUCTION 

Thousands of sun-grazing comets have been de- 
tected since the advent of space-based solar observa- 
tion. The most successful instrument for observing these 
has been the Lar ge Angle and Spectro metric Coron- 
agraph (LASCO; ferueckner ~a l. 1995) on the Solar 
and Heliospheric Observatory (SOHO; iDomingo et al.l 
1995), which has observed more t han 2000 such comets 
approaching the S un (see, e.g., iBiesecker et al.l [20021 : 
iKnight et al.l 12010). For a complete catalogue of sun- 
grazing comets, see the US Naval Research Labora- 
tory Sungrazing Comets websiteQ Very few sun-grazing 
comets have been seen to survive their close passage 
to the solar p hotosphere and emerge post perihelion 
(|Marsdenl (2005). none of which have been observed by 
LASCO until 201 10 

Despite the lack of direct observation of the destruc- 
tion of a sun-grazing comet, the literature is not lack- 
ing in predictions of the results of such an event. It is 
not possible to give a full literature review he re, but we 
highli ght so me paper s of pa rticular relevance. IWeissmanl 
(1983) and Sckanina (1984) discussed the physical pro- 
cesses of import to the destruction of comets during their 
perihelion passage. A more complete model det ailing the 
erosio n of sun-grazing comets was presente d by Sckanina 
( 2003). And, most recently, iBrown et ail (|2011[ ) have 
provided analytical models describing the destruction 
mechanism of sun-grazing and -impacting comets. This 
latter work indicates that the dominant mass-loss mech- 
anism varies between sublimation, ablation, and explo- 
sion depending on the cometary mass and perihelion dis- 
tance. One particularly interesting prediction is that a 
sun-impacting comet large enough to reach the chromo- 
sphere would result in solar flare-like energy release. 

Further from the Sun, it has been known for around 15 
years that comets emit at X-ray and EUV wavelengths. 
In 1996, comet C/1996 B2 (H yakutake) was ob served 
by the ROSAT X-ray telescope (jLisse et al.lll996ft . Sub- 

1 http://sungrazer.nrl.navy.mil/ 
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X-rays 



all 



comets withi n 3 AU were found to emit 
19971 ) . The discovery of this emis- 



(|Dennerl et al.l 
sion was initially surprising, given the low temperature 
of the cometary atmosphere, but was explained as aris- 
ing from charge exchange between solar wind ions and 
neutr al cometary species (|Cravensl Il997t iKrasnopolskvl 
1 1997T ) - However, until 2011, there had been no observa- 
tions of a sun-grazing comet close to perihelion, let alone 
at such wavelengths. 

The launch of th e Solar Dynamics Observatory (SDO; 
iPesnell et aLl[2012h in 2010 provides a unique capability 
to observe sun-grazi ng comets. The Atm ospheric Imag- 
ing Assembly (AIA; iLemen "et~aT1[20Tl images the en- 
tire Earth- facing solar corona at a cadence of 12 s with 
~ 1.2 arcsec resolution in 7 EUV channels, with pass- 
bands centered at 94, 131, 171, 193, 211, 304, and 335 A. 
Should a sun-grazing comet be sufficiently bright in the 
EUV as it flies within the field-of-view of AIA, the high 
spatial and temporal resolution of AIA makes it the ideal 
instrument to observe such an event. 

The first detection of a comet by SDO was on 2011 
July 5-6, when comet C/2011 N3 (SOHO) passed across 
the disk of the Sun (jSchriiver et al.ll2012h . It was de- 
tected in several of the AIA passbands before disappear- 
ing as it evaporated in the solar atmosphere. A second 
comet, C/2011 W3 (Lovejoy), was observed on 2011 De- 
cember 15-16, this time passing behind the east limb of 
the solar disk and emerging on the west limb. Again, 
several of the AIA passbands detected EUV emission. 
As we will show below, cometary neutral species cannot 
survive very long at such close proximity to the solar ra- 
diation field. Thus, the explanation for EUV emission in 
sun-grazing comets must differ from the charge-exchange 
model used to explain the X-ray emission detected in the 
heliosphere. This paper offers such an explanation. 

The remainder of the paper is organized as follows. 
In Section [2] we describe the AIA observations of the 
two comets. In Section [3] we describe the model of the 
comet-corona interaction used to explain the EUV emis- 
sion, including our simplified description of the cometary 
topology and the physical parameters of the cometary 
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and solar atmospheres used in our calculations. Section^] 
outlines the rate of evolution of the cometary material as 
it is sublimated from the comet body and subsequently 
dissociated and ionized in the corona. In Section [5] we 
show what contribution the cometary material has to 
EUV emission in the AIA channels. Finally, in Section |6] 
we give our conclusions and offer some suggestions on 
future work to improve our emission model. 

2. OBSERVATIONS 
2.1. Comet C/2011 N3 (SOHO) 

Comet C/2011 N3 (SOHO) was first detected by the 
LASCO coronagraph in white light, approaching the west 
limb of the Sun. On 2011 July 5-6 it pa ssed across the 
disk of the Sun and was detected by AIA (Schrijvcr et al. 
12012ft . The first detection was in AIA 171 A, off the so- 
lar limb, at 2011 July 5 23:46 UT, and was tracked until 
2011 July 6 00:05:50 UT. Emission was also detected in 
the 131, 193, 211, and 335 A passbands. The observa- 
tions suggest that the comet may have fragmented before 
evaporating in the solar atmosphere. 

Of the 7 EUV channels of AIA , 5 show a response to 
the comet (|Schriiver et all 12012ft . There was no signif- 
icant detection in the 94 or 304 A channels; the 94 A 
channel has a low signal-to-noise ratio in these observa- 
tions and the 304 A channel is dominated by the bright 
emission from lower in the solar atmosphere. The en- 
hancement over the background emission is greatest in 
the 131 and 171 A passbands, where the cometary emis- 
sion is ~ 10% above the background. However, the comet 
brightness was determined by summing over a box of 
30 x 15 AIA pixels — a significantly larger area than the 
emitting region — so the true intensity enhancement due 
to the comet is greater. 

2.2. Comet C/2011 W3 (Lovejoy) 

A second comet, C/2011 W3 (Lovejoy), was observed 
by SDO on 2011 December 15. This passage of another 
Kreutz-family comet 5 months later resulted in perihe- 
lion being on the far side of the Sun from the perspective 
of SDO. Ground observations several days prior to per- 
ihelion gave an accurate estimate of the cometary orbit 
(|Greenll2011ft . allowing SDO to repoint 1250 arcsec east 
of solar center, giving extended observations of the comet 
as it approached the eastern limb. Comet C/2011 W3 
(Lovejoy) was appreciably more massive than C/2011 N3 
(SOHO), allowing it to surviv e perihelion and emerg e on 
the other side of the Sun (see IGundlach et~aT1 (2012) for 
a discussion on the survival) . SDO had repointed to Sun 
center in time to observe the emer gence. An in-depth 
analys is of the orbit can be found in Sekanina fc ChodasI 
(2012). We make particular note of the absence of ob- 
servable cometary dust at the time of the AIA obser- 
vations as a result of any dust parti cles being rapidly 
sublimated ()Sekanina fc Chod"asll2012D . 

Despite perihelion passage being obscured by the solar 
disk, C/2011 W3 (Lovejoy) resulted in, perhaps, more 
interesting observations than C/2011 N3 (SOHO). The 
absence of contamination from background emission from 
the solar disk in the off-limb observations gives a clear 
picture of sublimated cometary material following the so- 
lar magnetic field. In contrast to C/2011 N3 (SOHO), 



C/2011 W3 (Lovejoy) produced significant signals in the 
304 A channel. For the channels that also detected 
C/2011 N3 (SOHO), the brightness contrast is similar 
in the C/2011 W3 (Lovejoy) observations. 

Fig. Q] shows the comet ingress in the 171, 131, and 
304 A wavelength channels of AIA. Two images at each 
wavelength are shown, separated by 1 minute, and nearly 
simultaneous across the three wavelength channels. Pan- 
els (a) and (b) show the 171 A emission, where evidence 
of the cometary material following the magnetic field 
is most clearly visible. The material is seen to form 
clear striations, roughly perpendicular to the direction 
of the comet's motion. These striations are long-lived; 
the brightest was first detected in the 171 A channel at 
2011 Dec 15 23:57:47 UT and seen to persist until 2011 
Dec 16 00:16:23 UT, almost 19 minutes, and 8 minutes 
after the comet went behind the solar limb. 

Panels (c) and (d) of fig.[T]show the 131 A observations 
at close to the same time as those of (a) and (b). The 
signal-to- noise is smaller than in the 171 A passband, but 
the cometary emission is spatially similar. This is not the 
case for 304 A emission, shown in panels (e) and (f). The 
emission shown in panel (e) is more spatially confined 
than that of (a) and (c). Comparing the emission at the 
later observation time shown in panel (f ) with (b) and (d) 
makes this even more apparent. For ease of comparison 
across the images, we have highlighted with an arrow a 
point where the emission from the comet tail appears to 
kink due to the magnetic field of the corona. Sunward of 
this point, cometary material appears to move up and left 
as it follows the magnetic field. In contrast, the emitting 
material below the arrow moves down and right. We 
also note that there is no detectable signal in the 304 A 
channel at this point at the time of these images. 

The emission in the 171 and 131 A channels is far 
longer-lived than that of 304 A. Under normal Quiet Sun 
(QS) conditions, these 3 channels image plasma at logT 
(K) of 5.9, 5.6, and 4.7 respectively. The similarity in the 
emission from the 131 and 171 A channels thus indicate a 
different ionization/excitation mechanism from the QS. 
We suggest a solution to this problem in Section [5j 

The differences in the comets' detections clearly indi- 
cate that characteristics of the comet and viewing geom- 
etry significantly influence the intensity of EUV emission 
caused by their interaction with the solar atmosphere. It 
is our aim in this work to propose a mechanism whereby 
emission is possible; characterizing the precise nature of 
the emission depends on an accurate knowledge of the 
cometary and coronal parameters. 

3. EMISSION MODEL 

As the comet approaches the Sun, material is subli- 
mated from the comet surface. The mass loss rate and to- 
tal mass lost of C/2011 N3 (SOHO) during the visibility 
in the AIA images is estimated to be (0. 01 - 1) x 10 8 g/s 
and ( 0.06 - 6) x 10 10 g, respectively, by iSchrijver et al.l 
(|2012f) . 

We have constructed a simplified model that describes 
the evolution of the cometary material as it interacts with 
the solar atmosphere. We represent the comet coma as a 
cylinder along the direction of travel. A schematic rep- 
resentation is shown in Fig. [2] The body of the comet is 
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Figure 1. Emission from comet C/2011 W3 (Lovejoy) pre-perihelion. The color scale has been reversed, with darker colors indicating 
brighter emission. The comet is traveling up and to the right, towards the solar disk, which is to the extreme right of the images. All panels 
show the same location. Axes units are given in AIA pixels with (0,0) being the bottom left of the full 4096 X 4096 AIA image. Panels (a) 
and (b) show emission from the 171 A channel of AIA at 00:00:12 UT and 00:01:11 UT, respectively; panels (c) and (d) show the 131 A 
channel at 00:00:10 UT and 00:01:10 UT, respectively; and panels (e) and (f) show the 304 A channel at 00:00:08 UT and 00:01:08 UT, 
respectively. The arrow is at the same location in each image and indicates the kink in the tail described in the text. 



represented by the red circle at the center with radius r c . 
We use r c = 50 m, the upper limit of the nucleus radius 
given bv lSchriiver et all (|2012f) . The material sublimated 
from the surface then expands into the corona and forms 
the cylindrical geometry shown, with axis along the path 
of the comet's motion. The atoms that are formed in the 
coma become more highly ionized as they expand into 
the corona. We represent the space occupied by these 
ions as concentric shells (shown in blue in Figure [2]) with 
inner radius r\ and outer radius r<i. These radii are de- 
termined by the dissociation and ionization rates and the 
outgassing velocity. 

The biggest failing of this model is that we do not take 
account of the magnetic field. Rather than a simple ra- 
dial expansion, once the cometary material is ionized, it 
will be affected by the coronal magnetic field. As the level 
of ionization increases, the Larmor radius of the ions will 
decrease and the more highly charged species will form 
tighter spirals around the field lines. This is somewhat 
in opposition to the model presented here. The obser- 
vations of C/2011 W3 (Lovejoy) clearly show that the 
cometary material is influenced by the coronal magnetic 
field, but whether the passage of the comet has any sig- 



nificant bearing on the underlying field has yet to be de- 
termined. Accounting for the magnetic field is not trivial 
and we leave such a study to a future work. However, 
the ions will move in a cylindrical pattern. 

The precise composition of the material sublimated 
from the comet's surface is unknown. For the purposes of 
this work, we adopt the elemental abundances of Co met 
lP/Halley (given in Table [J) from iDelsemmd (jl988H . H 
and O are the most abundant elements due largely to 
the presence of water ice in the comet and oxides being 
dominant in the asteroidal dust particles. Once it is de- 
posited within the corona, H is rapidly ionized and thus 
does not contribute in any significant degree to the EUV 
emission. The most important elements to consider are 
O and Fe; O has a high abundance in the comet and 
many of the AIA channels cover wavelength ranges with 
strong Fe emission lines. To be as complete as possi- 
ble, we have examined the 8 most abundant elements 
of Comet lP/Halley, i.e., those listed in Table [T] — H, C, 
N, O, Mg, Si, S, and Fe. We will show that, of these, 
only O and Fe cause any significant emission in the AIA 
passbands. 

There is some uncertainty in the outflow velocity from 
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Figure 2. Schematic of the cometary emission model. The figure on the left shows the cross section perpendicular to the direction of 
motion. The figure on the right shows the comet trajectory as seen by AIA. The red dot in the center indicates the comet body. The blue 
area indicates one of the shells described in the text where a given ion exists. 



Table 1 

Elemental Abundance of Comet 



Element Fractional Abundance (Number) 



H 0.484 

C 0.137 

N 0.023 

O 0.304 

Mg 0.011 

Si 0.016 

S 0.010 

Fe 0.011 



Note. — Elemental abundances in Comet lP/Halley 



JDelsemme!lT988T) . 



the surface; from the observations of C/2011 N3 (SOHO), 
we estimate the time taken for material to travel from the 
surface of the nucleus to the extent of the visible coma to 
be ~ 60 s. Assuming the material hows radially from the 
nucleus at a constant velocity gives an outflow velocity 
of w ut = 17 km s _1 . We explore the sensitivity of the 
model to this time scale by allowing it to vary from 10 s 
to 90 s. The emission predicted by the model described 
here is dependent on the total time the cometary mate- 
rial remains emitting, rather than on the precise outflow 
velocity, so any influence of the magnetic field on the di- 
rection and magnitude of v out will not greatly affect the 
model results if the emitting time remains unchanged. 
For comet C/2011 W3 (Lovejoy), however, emission is 
seen to persist for up to 20 minutes in certain AIA chan- 
nels. The nature of this emission is somewhat different 
from the simplistic model presented here. It is seen to 
follow the solar magnetic field but, after the first minute 
of emission, does not appear to expand into the corona, 
rather maintaining consistent spatial dimensions before 
gradually fading from view. For this reason, we have cho- 
sen to use the outflow velocities of v ou t = 17 km s _1 for 
the emission calculations, but we make some comments 
on the extended duration of the C/2011 W3 (Lovejoy) 
emission in Section [5] 

The comet periheli on parameter of q = 0.0052986 
(|Schriiver et all [20T1 for C/2011 N3 (SOHO) corre- 
sponds to a height above the solar surface of 97,200 km. 
Both comets were observed to pass through regions of low 
solar activity. Typical values for the QS electron temper- 
ature and density at this height are T c — 1.5 x 10 6 K and 
10 s cm -3 , respectively. We use these values through- 
out for the zone of interaction between the comet and 



corona. In the absence of the comet, emission measured 
by AIA comes predominantly from lower in the atmo- 
sphere. When comparing observations of the comet with 
background QS emission, we use an average QS differen- 
tial emission measure (DEM) from iVernazza fe Reevei 
(|1978f) and an electron density of 5 x 10 8 cm~ 3 . 

4. DISSOCIATION AND IONIZATION OF 
COMETARY MATERIAL 

4.1. Dissociation 

The subl i mated molecules are 50% water by mass 
(|Delsemmel I1988D . The destruction of these water 
molecules is mostly by photodissociation: 

H 2 + lw^OH + H + 3.42 eV (1) 

^H 2 + 0( 1 L>) + 3.84eV (2) 

-^H + H + O + 0.7 eV (3) 

-^H 2 0+ +e + 12.4 eV. (4) 

The first of these reactions is the most likely, with a 
branching ratio of 86%. The reaction r ate of this process 
at 1 AU is given bv lCombi et~aLl ()2004D as 1.04 x 10" 5 s^ 1 
for the nonflaring Sun at m edium activity. The rates 
given bv lCombi et a l. (2004) are in broad agree ment with 
those of previous work bv lHuebner et al.l 1)19921 ) . If we as- 
sume this rate scales as the inverse square of the distance 
to the radiation source, then the photodissociation rate 
of the comet at the solar surface is~4.8xl0 _1 s _1 . This 
translates to a lifetime of ~ 2.08 s for a water molecule at 
the solar surface. At small heliocentric distances, ther- 
mally excited states in the molecule will increase the pho- 
todissociation rate such that the r~ 2 scaling is not en- 
tirely accurate. Our estimate of the lifetime is thus an 
upper limit (i.e., a lower limit of the rate constant). 

The hydroxyl radical resulting from the dissociation of 
water will also photodissociate: 

OH + h^^O( 3 P)+H + 1.27eV (5) 

-^OCD) +H + 7.90 eV (6) 

-^OCS) + H + 9.80 eV (7) 

-^OH++c + 19.1eV. (8) 

To find the total rate for OH dissociating to O and H, we 
take the sum of the first three of the se reactions. The to - 
tal rate constant at 1 AU is given bv lCombi et "a l. (2004) 
as 2.5 x 10 -5 s _1 . The resulting rate coefficient at the so- 
lar surface is then 1.2 s , corresponding to a lifetime of 
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0.86 s. The mean time for a water molecule to dissociate 
to O and H can then be estimated as 2.9 s. 

In addition to water, various other molecules are subli- 
mated from the surface. Under typical solar conditions, 
most of the AIA channels are dominated by emission 
from Fe ions so we are interested in any Fe sublimated 
from the comet that may contribute to the emission de- 
tected by AIA. Fe I emissio n has been ob served in the 
spectra of comet Ikeya-Seki ()Prestonl 119671 ) . most likely 
produced by ablation or vaporization of refractory grains. 
At the end stages of the asteroidal evaporation, the most 
important molecule containing Fe is ferrous oxide (FeO). 
The photodissociati on cross section for FeO at 252 nm 
is 1.2 x 10" 18 cm 2 (|Chestakov et al.ll2005D . The solar 
flux at 1 AU is ~ 5 x 10 12 photons cm~ 2 s _1 at 252 nm, 
so the rate at 1 AU is at least ~ 6 x 10~ 6 s" 1 . At the 
solar surface this is 0.28 s" 1 , or a lifetime of 3.6 s, i.e., 
comparable to that of water. Other molecules, such as 
MgO, will also rapidly photodissociate, adding another 
O atom to the cometary debris. However, the metallic 
ion so produced does not emit in the AIA passbands. 

4.2. Ionization 

The neutral atoms formed by dissociation are ionized 
through charge exchange with coronal protons, through 
impact with coronal electrons and protons, and by pho- 
toionization. We consider each of these processes in turn. 

4.2.1. Charge exchange 

Electron capture by hydrogen ions due to collisions 
with neutral oxygen, 

H+ + 0( 3 P) ^H + 0+( 4 S'), (9) 

is a near res onant reaction . Cro ss sections have been 
calculated by iStancil et al.l (|1999f ) for collision energies 
between 0.1 meV/u and 10 MeV/u. For a collision en- 
ergy of 2 keV (corresponding to the energy of H + in the 
collision with the comet material) the rate coefficient for 
the reaction is 5.52 x 10~ 8 . For a coronal proton den- 
sity of 10 s cm -3 , the lifetime of cometary O atoms is thus 
0.18 s. The reverse reaction is considerably less likely due 
to the negligible abundance of neutral H in the corona. 

The resonant charge exchange process between H and 
protons, 

H+ +H(ls) -> H(ls) +H+, (10) 

has been studied bv lBates fe Dalgarnol (|1953l ). They pro- 
vide tabulated cross sections down to an impact energy 
of 2.5 keV. The cross section has a weak energy depen- 
dence and no other resonances at these energies so we 
can extrapolate below this energy to find a cross section 
of ~ 6 x 10" 15 cm 2 at an energy of 2 keV. For a coro- 
nal proton density of 10 8 cm -3 and a collision energy of 
2 keV, this translates to a mean lifetime of only 0.027 s 
before cometary hydrogen atoms are ionized by charge 
transfer. As above, the inverse process is much less effi- 
cient because of the scarcity of neutral H in the corona. 

During the initial expansion the density of H and O 
exceeds 10 8 cm" 3 . This means the time to become ion- 
ized is increased until the comet material has expanded 
to a radius that includes sufficient protons to ionize all 
H and O. Until then, the situation is more complicated. 



4.2.2. Electron-impact ionization 

Following the initial ionization of the neutral species, 
subsequent ionization is dominated by impact with free 
electrons of the corona. We estimate the timescales of 
these reactions under the assumption that the electron 
velocity is large in comparison to that of the cometary 
ions. For a Maxwellian distribution of free electrons with 
temperature 1.5 x 10 6 K, the mean electron velocity is, 



v e = W— ^ « 7.6 x 10 3 kms -1 , (11) 
V nm c 

so the free electrons in the corona are over an order 
of magnitude faster than the cometary oxygen ions. 
For comparison, coronal protons and oxygen ions have 
mean velocities of 1.8 x 10 2 and 44 km s _1 , respec- 
tively. A proper analysis of the electron-impact ioniza- 
tion rate would consider an anisotropic electron veloc- 
ity distribution in the rest frame of the ion. However, 
given the disparity in the velocity of the interacting elec- 
trons and ions, it is safe to proceed under the assump- 
tion of an isotropic Maxwellian electron distribution of 
T c = 1.5 x 10 6 K in the rest frame of the target ion. 

To trace the ionization level of cometary oxygen as it 
passes through the corona, we consider the reactions: 

O q+ + c -> (9+1)+ + e + e, (12) 

where q is t he charge of the ion. We use the rate 
coefficients of Der^ (|2007f) . an electron temperature of 
1.5 x 10 6 K, and a density of 10 8 cm to derive reac- 
tion rates for this process. The results are summarized 
in Tabled 

From these results we can estimate which ionization 
stages of O will be sufficiently long-lived to contribute 
to EUV emission. If the plasma were allowed to re- 
lax to ionization equilibrium then the highest ioniza- 
tion stage reached u nder these conditions would be 6+ 
(jBrvans et al.l 120091 ) . We thus expect to see emission 
from ionization stages up to and including q = 6. Given 
the ionization rates of the cometary material, this ion- 
ization state will be reached only after sufficient time has 
el apsed . 

IDerd (|2007j) also provide rate coefficients for the 
electron-impact ionization of Fe. We have given the rates 
for the first 12 ionization stages at an electron tempera- 
ture of 1.5 x 10 6 K and a density of 10 8 cm" 3 in Tabled 
Analyzing these timescales shows that the highest ion- 
ization stage reached is Fe 9+ for v out = 17 km s _1 . 

4.2.3. Photoionization 

The photoionization rate of the dissociated atoms and 
ions is dependent on the solar EUV flux. The comet 
is exposed to a wide range of radiances as it travels 
close to the solar surface. For this reason, we choose 
to use the average EUV spectral irradiances at 1 AU 
from the Thermosphere Ionosphere Mesosphere Energet- 
ics Dynamic s (TIMED) spacecr aft's Solar EUV Exper- 
iment (SEE; iWoods et al.l 120051 ). Values between 2002 
and 2011.5 were averaged to produce the irradiances used 
here. We use the cross section s from lVerner fe Yakovlevl 
(1995) and lVerner et al.l (|1996l ) to calculate the photoion- 
ization rates at 1 AU and have calculated the equivalent 
rate at the solar surface by scaling by r 2 . The values are 
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Table 2 

Ionization Rates for Selected Ions (s — 1 ) 



Ionization Level 


CX 


EII 


Photoionization 


PII 


Total 


H -> H+ 


37.6 


2.9 


4.3 X 10~ 3 


~ 


40.5 


O -y 0+ 


5.5 


8.3 


1.5 X 10~ 2 


4.2 x 10 _1 


14.2 


o+ -> o 2+ 


— 


2.8 


6.5 x 10~ 3 


3.1 x 10~ 4 


2.8 


O + —y O + 




1.2 


1.4 x 10 -3 


4.2 x 10~ 5 


1.2 


— > o + 




4.1 X 10 _1 


5.1 x 10 -5 


8.0 x 10 -6 


4.1 x 10 _1 


O a-r _^ Q oi- 




1.1 x 10 _1 




1.5 x 10 -6 


1.1 x 10" 1 


O aT — y O t 




3.4 x 10 -2 




3.5 x 10~ 7 


3.4 x 10~ 2 


Fe — ► Fe"' - 




20.6 


8.4 X 10 


1.7 


22.3 


Fe + — > Fe 2+ 




6.3 




3.6 x 10~ 3 


6.3 


Fe 2 + -j- Fe 3+ 




6.6 




5.4 x 10~ 4 


6.6 


Fe 3 + ->• Fe 4 + 




3.8 




8.2 x 10~ 5 


3.8 


Fe 4 + -> Fe 5 + 




1.6 




1.9 x 10~ 5 


1.6 


Fe 5 + -y Fe 6+ 




8.4 x 1CT 1 




5.5 x 10" 6 


8.4 x lO" 1 


Fe 6 + -> Fe 7 + 




4.6 x 1CT 1 




1.9 x 10" 6 


4.6 x lO" 1 


Fe 7 + -> Fe 8 + 




2.0 x 10" 1 




6.8 x 10~ 7 


2.0 x lO" 1 


Fe 8 + -> Fe 9 + 




5.4 x 10" 2 




3.0 x 10" 7 


5.4 x 10~ 2 


Fe 9 + -> Fe 10 + 




3.0 x 10~ 2 




1.7 x 10- 7 


3.0 x 10~ 2 


Fc 10 + -> Fc 11+ 




1.5 x 10~ 2 




9.5 x 10~ 8 


1.5 x 10" 2 


Fe n + -> Fe 12 + 




7.5 x 10" 3 




4.8 x 10" 8 


7.5 x 10~ 3 


Fe 12 + -> Fe 13 + 




4.0 x 10~ 3 




2.6 x 10~ 8 


4.0 x 10" 3 



Note. — Ionization rates of various processes at the solar surface. Charge exchange (CX), electron-impact (EII) and proton- 
impact ionization (PII) rates are for an electron and proton density of 10 8 cm -3 . Photoionization was not considered for states 
4+ and higher, or Fe + and higher. 



given in Table [5J In comparison to resonant charge ex- 
change and electron-impact ionization, photoionization 
has a negligible impact on the ionization rate. A flare 
could have larger photoionization consequences, but none 
were observed during the perihelion passage of either of 
the comets discussed here. For this reason we have only 
listed the photoionization rate for the low charge states 
of O and Fe in Table H 

Photoionization of neutral hydrogen al so plays a rol e 
in the production of protons. According to lKellerl (|1976f ). 
the calculated lifetime of cometary H atoms due to pho- 
toionization from the solar radiation field at a distance 
of 1 AU is 1.4 x 10 7 s. This value is comparable to that 
from the TIMED/SEE irradiances, where the average 
photoionization lifetime over 2002-2011.5 was calculated 
as 1.1 x 10 7 s at 1 AU. The resulting lifetime against 
photoionization at the solar surface is 2.4 x 10 2 s. 

4.2.4. Ionization by protons 

The relative velocity of the cometary ions to the pro- 
tons in the corona is ~ 600 km s _1 . To calculate the rate 
of ionization due to impact wit h energetic proto ns we use 
the classical approximation of iGrvzin ski (1965). Results 
are shown in Table [2] for O and Fe. They are negligible 
compared to other ionization processes. 

Cometary protons, resulting from sublimated water 
molecules that have been dissociated and subsequently 
ionized, will initially be traveling at the velocity of the 
comet relative to the corona. This velocity of 600 km s _1 
corresponds to a proton beam of 2 keV. Can these fast 
protons cause ionization and excitation of the coronal 
plasma they are passing through? This is the same reac- 
tion as described above where fast cometary ions collide 
with slower coronal protons. Using the same approach, 
we can show that the proton impact does not have a 
significant effect on the ionization balance of the plasma. 



For it to be important, the proton-impact rate has to 
be comparable to the electron-impact ionization rate and 
the recombination rate. Taking ionization of O as an 
example, the rate coefficient for proton impact is ~ 5 
orders of magnitude smaller than that for electron impact 
for the ions present in the corona (see Table [2]). For the 
conditions outlined in Scction[3]and the dissociation time 
for water given in Section I4.1[ the highest proton density 
we can expect in the comet is ~ 10 11 cm -3 , only around 
3 orders of magnitude greater than the coronal density. 
The ionization rate by proton impact is thus 2 orders of 
magnitude smaller than that by electron impact. 

5. RESULTS 

Using the physical conditions outlined in Section [3] and 
the ionization rates of Section |4~21 we have calculated the 
line emission from all elements l isted in Table [p i n eac h 
AIA bandpass using CHIANTI (jDere et al.lll997l [20091) . 
Collisional redistribution among the excited states of the 
ion is several orders of magnitude faster than ionization, 
so we have assumed detailed balance applies to the dis- 
tribution of those energy levels. The results are shown in 
Figures [3H9] for an outflow time of 60 s. The intensities 
have been folded with the effective areas o f the respective 
AIA filters (taken from SolarSoft; see also lBoerner et all 
120121 ) and compared to the average QS emission. All in- 
tensities are integrated over the emitting area described 
in Section [3] and outlined in Fig. [2] and compared with 
QS intensities over the same area. Of the elements con- 
sidered (H, O, C, N, Mg, Si, S, and Fe), only the ions 
of O and Fe contribute any significant emissions in the 
wavelength ranges of the AIA instrument. 

We note here that the He n lines that dominate the 
QS emission in the 304 A channel are not properly mod- 
eled by CHIANTI, which underestimates the observed in- 
tensities by a significant factor (see, e.g.. lAndretta et all 
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170 175 
Wavelength (A) 

Figure 3. Emission from cometary ions in the bandpass of AIA 
171 A compared to QS. Fe emission from the comet is shown in 
red and O emission from the comet in blue. Black lines indicate 
a typical spectrum due to the QS. The dotted line is the effective 
area of the AIA passband. The colored lines are the spectrum from 
cometary O and Fe ions with outgassing velocity of 17 km s — 1 . The 
strongest emission lines have been labeled. 





Figure 6. As Fig. [3]but for 211 A. 



€ 0.6 
< 



S 0.4 



0.2 



He ! 



O III 



O IV 



220 



300 

Wavelength (A) 



128 130 132 

Wavelength (A) 



Figure 7. As Fig. \3\ but for 304 A. 



Figure 4. As Fig. [3] but for 131 
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Figure 5. As Fig. [3] but for 193 . 



Wavelength (A) 
Figure 8. As Fig. |3] but for 335 A. 




Figure 9. As Fig. [3] but for 94 A. 

2003). For this reason, we have multiplied the back- 
ground QS emission from these lines b y a factor of 2 0.0 
to bring them in line with observation ( Warren 2005). 

For each AIA bandpass we have calculated the frac- 
tional increase in radiance it would detect due to the 
cometary emission mechanism. We have integrated the 
emission over the wavelength range of each bandpass and 
compared the total intensity of the QS plus comet emis- 
sion to QS emission alone. We have calculated the de- 
pendence of this fractional increase on the outflow ve- 
locity and show the results in Fig. [TO] The time for the 
cometary material to flow from the nucleus to the edge 
of the emitting region is allowed the range 10-90 s (an 
outflow velocity of 11-100 km s _1 ). The solid lines show 
the contribution from O and the dashed lines show the 
contribution from Fe. The outflow velocity has an im- 
pact on the relative contribution from O and Fe ions to 
the emission spectra, most notably in the 131 and 171 A 
passbands, but does not alter the total fraction of emis- 
sion from cometary ions relative to the background to a 
significant degree. 

Our model predicts emission in the 171 A bandpass 
from O v, O VI, Fc ix, and a smaller contribution from 
Fe x. Under normal QS conditions this bandpass is 
dominated by Fe IX emission. The only other bandpass 
with significant emission from both O and Fe due to the 
cometary material is the 131 A channel. Here, contribu- 
tions to the emission come from O vi and Fe viii. These 
same Fe VIII lines are those normally strongest in this 
channel under QS conditions. It is significant that the 
171 and 131 A channels are the only two that detect 
strong O VI emission lines. These are the highest ioniza- 
tion stages of O detectable in any of the AIA passbands, 
taking longer to form than the lower stages, and thus 
explaining the time delay in the observations from these 
channels (see Fig. [1]). We expand on this argument in 
Section [5] 

In the 193, 211, 304, and 335 A passbands, O emis- 
sion dominates over Fe. A single O v line dominates the 
emission spectrum in the 193 A passband. The 211 A 
passband has contributions from a number of O iv and 
O v lines. O in and O IV are the strongest emitters in 
the 304 A passband, and O in dominates the emission in 



Figure 10. The intensity of O and Fe emission as a function of 
outflow velocity as detected by each AIA bandpass is given as a 
fraction of the total emission from the comet plus the average QS 
emission. The results are given for the conditions given in Sec- 
tion [3] Solid lines show the fractional contribution from cometary 
O ions and dashed lines from cometary Fe ions. 

the 335 A passband. We note a progression of emission 
from O v, through O iv, to O in on advancing through 
the 193, 211, 304, and 335 A passbands. We return to 
this point in Section [5] when discussing a possible time 
lag in the observations from these channels. 

The 94 A passband has no O emission lines in its wave- 
length range. Here, emission from the comet material is 
from Fe viii and Fe x and acts to increase the same emis- 
sion lines that are strongest in the quiescent QS. 

6. DISCUSSION AND CONCLUSIONS 

We have shown that emission from the O and Fe atoms 
deposited in the corona by the passage of a comet is 
a viable mechanism to explain the increased radiance 
detected by AIA. The rapid injection of atoms to the 
corona followed by their ionization through successive 
ionization stages produces emission lines not normally 
present at coronal temperatures. The exact structure of 
the resulting spectra and their magnitude relative to the 
background corona are dependent on both the coronal 
and cometary parameters. Comparing our calculations 
with the AIA observations may provide a means of deter- 
mining these parameters more accurately than has been 
possible heretofore. 

The difference in the emission detected across the AIA 
channels for C/2011 W3 (Lovejoy) can be explained by 
our model if one considers the O emission. There is some 
observational evidence of the 304 and 335 A channels 
showing emission only close to the comet nucleus but that 
does not persist. By contrast, the 131 and 171 A channels 
show more spatially extended emission that persists for 
up to 20 minutes. The 304 A channel has contributions 
from O in and O IV while 335 A has strong contributions 
from O iv. According to our model, O will emit in these 
ionization stages for several seconds before being further 
ionized. Both 131 and 171 A channels have significant 
contributions from O VI. This ion has a higher formation 
temperature, close to the ambient coronal temperature, 
and so remains present in the plasma for longer than the 
lower ionization stages. The observed fading must come 



9 



from the diminution in the number density of oxygen 
nuclei in the field-of-view. 

By a similar argument, the model presented here pre- 
dicts that there should be a time lag in the emission 
from the AIA channels. This should follow the ioniza- 
tion stages of O that contribute to the emission in each 
wavelength channel, i.e., emission should first be detected 
in 335, followed by 304, then 211 and 193, and finally in 
131 and 171 A. Lightcurves of C/2011 W3 (Lovejoy) are 
being measured and the time lag between the channels 
is expected to be only a few seconds (see Tabled]) so it 
remains to be seen whether observations spaced by 12 
seconds can resolve any such time delay. 

Estimates of several parameters influence the predicted 
emission to varying degrees. The assumed density of 
the corona through which the comet passes is partic- 
ularly important. A larger density than that assumed 
here (10 8 cm' 3 ) would result in faster ionization and 
hence increased emission from higher ionization stages of 
O and Fe. The density assumed for the plasma respon- 
sible for the QS emission does not affect the cometary 
emission but will change the spectrum of the background 
corona that interferes with the detection of the comet. 
Finally, the density of the outgassing molecules will im- 
pact the relative intensity of the cometary and QS spec- 
tra. Within the range of reasonably expected conditions 
(around l-2xl0 6 K), the temperature assumed for the 
corona does not significantly alter the calculated spectra. 

Our model assumes the cometary species move radially 
away from the nucleus of the comet after sublimating. 
This assumption neglects the influence of the magnetic 
field on the ions. The primary effect of the gyrorotation 
will be an increase in the effective pathlength of the ions 
as they spiral along the field lines. A proper analysis 
would require the direction and magnitude of the mag- 
netic field, and is beyond the scope of this work, but is 
worth further study. 

Intensity profiles from AIA channels also show indi- 
cations of absorption by the denser inner regions of the 
comet's coma when v iewed against the disk of the Sun 
(|Schriiver et al.l [2012T ). The model presented here may 
be able to account for this. Near the comet surface, ions 
have yet to be ionized to the stages responsible for the 
emission seen later in the AIA observations. It is possible 
that the lower ionization stages cause absorption in this 
region and the relative brightness may track the O I pho- 
toionization curve. We intend to explore this possibility 
in future work. 

Research was supported by the Solar Dynamics Obser- 
vatory. 
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